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Abstract 

We study the ground state of a classical X-Y model with p > 3-fold spin 
anisotropy D in a uniform external field, H. An interface is introduced into 
the system by a suitable choice of boundary conditions. For large D, as 
H — > 0, we prove using an expansion in D~ l that the interface unbinds from 
the surface through an infinite series of layering transitions. Numerical work 
shows that the transitions end in a sequence of critical end points. 

PACS numbers: 0550; 6470; 7540D 
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When an interface unbinds from a surface it can do so through a first or second order transi- 
tion or, on a lattice, via a sequence of first order layering transitions. In general the nature of 
the unbinding, or wetting, transition depends on the details of the microscopic interactions 
and external parameters such as the temperature and applied magnetic field [fj],@. 

In a discrete spin model at zero temperature long-range interactions are needed to sta- 
bilise series of first order layering transitions. However, if the interactions are short range, an 
infinite sequence of such transitions can occur at finite temperatures ||. This results from 
competition between the localising effect of the binding potential and the entropic terms in 
the free energy which favour a delocalised interface. 

In this paper we aim to discuss the role of a hithertofore unexplored parameter on the un- 
binding transition: the spin anisotropy. We shall show that, as discrete spins soften, layering 
transitions can be stabilised in simple, short-range clock models, even at zero temperature. 
An expansion in inverse spin anisotropy allows us to prove that an infinite sequence of 
layering phase transitions exist. Moreover, because the interesting features occur at zero 
temperature it is possible to follow the phase diagram numerically for all values of the spin 
anisotropy. In particular we are able to demonstrate how the boundaries between the dif- 
ferent interface phases end in critical end points and to pinpoint these with considerable 
precision. 

We consider the classical X-Y model in a magnetic field H with p-fold anisotropy D. 
The model is defined by the Hamiltonian 

N 

H = J2 {-Jcoa{9i-i ~ Oi) - H{cos6i - 1) - D (cos p6i - l)/p 2 } (1) 
i=i 

where i labels the spins on a one-dimensional lattice and 9% can take values between and 
2tt. An interface is imposed on the system by choosing 9 = 2ir/p, and 9^ = and letting 
N — ► oo. If the interface lies to the right of the n th spin from the end of the chain (that is, 
for #o fixed, between i — n — 1 and i = n) the corresponding interface phase will be labelled 
<n>. 

For infinite D the Hamiltonian ([[]) describes a p-state clock model. For H > the 
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interface is bound to the surface in state (1). The point (H = , D = oo) is a multiphase 
point where the interface has the same energy whatever its position on the lattice. As the 
anisotropy is reduced from infinity it becomes, for p > 3, energetically favourable for the 
spins to relax from their clock positions. This results from competition between the field 
term, which favours an interface in position (1) and the exchange interaction which prefers 
to minimise the angle between the spins and hence favours a free interface. The result of 
this is that the degeneracy at the multiphase point is broken and the interface unbinds from 
the surface through an infinite series of first order layering transitions (1), (2), (3) . . . . 

For large D the existence of these transitions can be proven using an expansion in D^ 1 [|J . 
Writing 

Oi = 0° + §i (2) 

where 0° = 6i(D = oo), and keeping only terms quadratic in the angles and their differences 
the Hamiltonian ([I]) becomes 

oo 

U = Yl {~ Jc i - Hcos{9? - 1) + Jq{^_i -§i + Si/dY/2 - Js 2 i /2c i + 
i=i 

Hcos9?{§i + tan0?} 2 /2 -if sin 2 0?/2 cos0? + DO? (3) 

where 

* = sm(*?-i-0|). * = 008(0^-0°). (4) 

The equilibrium values of the 9{ are given by minimising the Hamiltonian (|3p . This leads 
to linear recursion equations 
- Jci{9i-i -9~i + s t /a} + Jc i+1 {9i - 9 l+1 + s i+1 /ci +1 } + Hco$9° i {9 i + tan0°) + D9 { = 0. 

(5) 

If the full Hamiltonian ([]]) is used non-linearities appear in the recursion relations (0). 
However, these do not affect the leading order terms needed for the subsequent calculations. 

Leading order corrections to the values of the spins can easily be read off from the 
recursion equations (H). For the phase (n) 
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k+2 = (J/Dfsm(2ir/p) + 0(1/ D 4 ) 
~6 n+ i = (J/D) 2 sm(2ir/p) + 0(1/ D 3 ) 
O n = (J/D) sin(27r/p) + 0(1/ D 2 ) 

B n -x = -J/Dsm(2n/p) + 0(H/D, 1/D 2 ) 
O n -2 = -(J/D) 2 sin(27r/p) + 0(H/D, 1/D 3 ) 

~9 1 = -(J/D)"- 1 sin(27r/p) + 0(H/D, 1/D n ) (6) 

The final result for the interface phase boundaries derived below will demonstrate that it is 
consistent to neglect terns 0(H/D) in ([|). Note that for p = 2, for D large, the spins do 
not lower their energy by canting. 

Using the harmonic approximation (|^) the energy differences E/ n \ —E/ n -i\ between neigh- 
bouring interface states can be calculated. Let (n — 1) have spins {a,} with «i the surface 
spin and (n) have spins with ft the surface spin. Then in both cases the interface lies 
between i = n — 1 and i = n and a° = /3°. Using this labelling and the recursion equations 
(|) some algebra leads to the result 

E {n) - = -H {cos(2tt/j9) - 1} + Hfe sin $/2 - J (3^/2 n>3 (7) 

E {2) -E {1) = -H {cos(2tt/p) - 1} + Hf3 1 smp°/2 - J cos(2ir/p)f3 1 a 2 /2 + Jsm(2n/p)(3 1 /2 (8) 

These formulae are exact for the quadratic Hamiltonian (|3|). Higher order terms in the full 
Hamiltonian (|]) appear as higher order corrections. 

Substituting in the values for the surface spins from (Q) gives, for n > 2, 

E {n) - £ (n _D = -H{cos(2tt/ P ) - 1) - J 2n ~ 2 sin 2 (2n/p)/(2D 2n - 3 ) + 0(1/ ' D 2n ~ 2 ) (9) 

It follows immediately from (^|) that the boundary between phases (n — 1) and (n) is given 
to leading order by 
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#<„-!):<») = J 2 "' 2 sin 2 (27r/j9) {2 (1 - cos(27r/j9)) D 2n ~ 3 } (10) 

indicating that the unbinding proceeds via an infinite series of phases (n) of widths 
0(1/D 2n - 3 ). 

These results were confirmed numerically for the case p = 3 by studying iteratively the 
equations which minimize the energy ([1]). The numerical approach allowed us to obtain 
the interface phase diagram for all values of D which is shown in Figure 1. The first order 
boundaries between the different interface phases end at a series of critical end points at 



(1): 


:(2) 


D* 
-^1,2 


= 1.1268 ±0.0003 


TT* 

#1,2 


= 0.2566 ±0.0009 


(2): 


:(3) 


D* 


= 0.9360 ± 0.0003 


TT* 

#2,3 


= 0.04357 ±0.00001 


(3): 


:<4> 


D* 

u 3£ 


= 0.7281 ±0.0005 


TT* 

#3,4 


= 0.01029 ± 0.00004 


(4): 


:(5) 


D* 


= 0.5931 ±0.0003 


TT* 

#4,5 


= 0.00295 ± 0.00002 



These were identified as the points where both the energy E and its partial derivative with 
respect to H become the same in the two phases. Assuming that D* n+1 and H* n+1 have 
a power law dependence on n, (D^,!!^) = (0,0) is consistent with the data. However, it 
was only possible to obtain results for n < 5, and so we cannot be confident that this is the 
true asymptotic behaviour. 

For D = the interface shape varies continuously from being a domain wall at the surface 
for H large to a uniform spiral for H = 0. 

To summarise, we have shown that the softening of discrete spins provides a mechanism, 
somewhat analogous to temperature, which can stabilise interface layering transitions. An 
expansion in inverse spin anisotropy was used to demonstrate that an infinite number of such 
transitions exist at large D. Because the transitions take place at zero temperature it was 
possible to obtain good numerical estimates of the end points of the first order transition 
lines. 
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The model we describe is widely used in the theoretical description of magnetism in the 
rare-earth metals and compounds |5]||. The results presented here may be of relevance to 
attempts to model mutilayers comprised of these compounds 0j. Moreover, we hope that 
the system will have theoretical applicability in that it represents one of the simplest models 
where expansion about a multiphase point is feasible. Aubry || has recently pointed out 
that such multiphase points (which he calls anti-integrable limits) also exist in electronic 
systems and arrays of non-linear coupled oscillators. It may be possible to use extensions of 
the technique presented here to treat these systems. 
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Figure Caption 

Fig 1: Phase diagram of the classical X-Y model in a magnetic field, H, with 3- fold spin 
anisotropy, D, and an imposed interface. There are an infinite number of interface layering 
transitions. The first order boundaries between them terminate in critical end points. 
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